ggplot(data=df,aes(x=type,y=VAF,group=sample,color=group))+
geom_line(size=1)+geom_point(aes(shape=type),size=3)+
ylim(0,1)+theme_light()+xlab("")+ylab("variant allele frequency")+
theme(axis.title.y=element_text(size=14),axis.text.x=element_text(face="bold",size=12))

ggplot(data=c,aes(x=type,y=VAF,group=group,color=group))+
geom_line(size=1)+geom_point(aes(shape=type),size=3)+
ylim(0,1)+theme_light()+xlab("")+ylab("variant allele frequency")+
theme(axis.title.y=element_text(size=14),axis.text.x=element_text(face="bold",size=12))

file=read.table("result_HMM_enrichment.csv",head=T,sep=",")
groupf=unique(file$con.file.name)
file1=file[file$con.file.name==groupf[1],]
file1$HMM.group=paste0("NO.",file1$HMM.group)

anno=read.table("group意义.txt",head=T,sep="\t")
anno=data.frame(HMM.group=anno$STATE.NO.,description=anno$DESCRIPTION,num=anno$NUM)
filee=merge(file1,anno,by="HMM.group")
filee$`log2(OR+1)`=log2(filee$OR+1)
library(ggplot2)
 filee$group="nonsignificant"
 filee[filee$p.value<0.05,]$group="significant"
 

	
	ggplot(data=filee,aes(x=`log2(OR+1)`,y=num,group=group,color=group))+geom_point(aes(shape=group),size=3)+
    scale_y_continuous(breaks = filee$num,labels = filee$description)+xlim(0,3)+theme_light()+ylab("Term")+geom_vline(xintercept = 1,color="black")+
    scale_color_manual(values = alpha(c('green','red'),0.8))+theme(panel.grid.minor = element_blank())